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It is well known that magnetic activity in late-type stars increases with increasing rotation rate. Using inversion techniques 
akin to medical imaging, the rotationally broadened profiles from such stars can be used to reconstruct 'Doppler images' 
of the distribution of cool, dark starspots on their stellar surfaces. Interacting binaries, however, contain some of the 
most rapidly rotating late- type stars known and thus provide important tests of stellar dynamo models. Furthermore, 
magnetic activity is thought to play a key role in their evolution, behaviour and accretion dynamics. Despite this, we 
know comparatively little about the magnetic activity and its influence on such binaries. In this review we summarise the 
concepts behind indirect imaging of these systems, and present movies of the starspot distributions on the cool stars in 
some interacting binaries. We conclude with a look at the future opportunities that such studies may provide. 
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1 Introduction 

Starspots, the most easily observed manifestation of mag- 
netic activity, are a ubiquitous feature of rotating late-type 
stars with outer convective zones. Indeed, naked-eye ob- 
servations of spots on our own Sun first took place lOOO's 
of years ago, and subsequent telescopic observations have 
been undertaken almost continuously for over 3 centuries. 
Despite this, our understanding of how magnetic fields are 
produced and evolve on the Sun and other stars, and how 
fundamental parameters affect magnetic field generation, is 
still critically lacking. Since Vogt & Penrod (1983) demon- 
strated that starspots could be seen as bumps traversing the 
Une profiles of HR 1 099, several indirect imaging techniques| 
have been used to map the magnetic topographies on mainly 
(but not exclusively) isolated rapid rotators. Since it is well 
known that magnetic activity increases with increasing stel- 
lar rotation rate, much of this work has focused on deter- 
mining how rotation impacts starspot emergence. 

Interacting binaries such as Cataclysmic Variables (CVs)| 
contain some of the most rapidly rotating late-type stars 
known. CVs typically consist of a late-type main-sequence 
star in orbit around a more massive white dwarf. The close 
proximity of both components causes the late-type (or donor)| 
star to fill its Roche-lobe and transfer mass to a white dwarf. 
With orbital periods as short as an hour, and tides forcing 
synchronous rotation of the donor, images of these systems 
provide excellent tests of stellar dynamo theories at ultra- 
fast rotation rates. In addition, such systems also allow the 
impact of tides on the stellar dynamo to be assessed. 
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Furthermore, many interacting binaries show variations 
in orbital periods, brightnesses, and outburst durations that 
have been attributed to starspots and activity cycles. Despite 
their obvious importance, both to understanding stellar dy- 
namos and understanding the observed behaviour of inter- 
acting binaries, little progress in this field has been made 
until recently. In this review, we will describe the astro- 
tomographic techniques required to image starspots on the 
cool stars in interacting binaries, before summarising some 
of the recent highlights in this field. Finally, we will con- 
clude with some future opportunities. 

2 The technique of Roche Tomography 

Unfortunately, interacting binaries are too small and too dis- 
tant to be resolved directly, even in the worlds largest tele- 
scopes. Rutten & Dhillon (1994), however, described a tech- 
nique, called Roche tomography, which uses phase-resolved 
spectra to reconstruct the line-intensity distribution on the 
donor stars. While analogous to Doppler imaging methods 
applied to single stars, Roche tomography differs in two 
important ways. First, the donor stars are tidally distorted 
into a Roche-lobe and are in synchronous rotation around 
the centre-of-mass of the binary. Second, the systems often 
show rapid variability due to the accretion. This means that 
one usually requires simultaneous photometry with which 
to slit-loss correct the spectra before continuum subtraction, 
whereas with Doppler imaging the spectra need not be slit- 
loss corrected and the continuum is normalised. 

Otherwise, the concept of Roche tomography is rela- 
tively straight forward. Given the binary parameters, the 
donor's surface can be defined to lie on the critical surface 
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of equal gravitational potential given by the Roche approxi- 
mation at the point where Roche-lobe overflow occurs. The 
surface can then be modelled as a series of tiles lying on this 
critical surface. Each tile is then assigned a copy of the local 
specific intensity profile (which may either be Gaussian or 
taken from a spectral-type standard) and convolved with the 
instrumental resolution. One can then assign an intensity to 
each tile which scales its contribution to the total profile. 
Each tile's contribution is then integrated over the visible 
surface of the star taking into account limb-darkening, fore- 
shortening and the radial-velocity of the elements. 

Roche tomography essentially carries out the reverse of 
the process just described. The contribution of each ele- 
ment is iteratively varied until a map is obtained which pre- 
dicts data consistent with the observed data. The goodness 
of fit is measured using the reduced statistic. Unfortu- 
nately, the reduced constraint alone is not sufficient to 
select a unique map as there are many maps that can fit the 
data equally well, and so we also adopt an additional reg- 
ularisation statistic. Following Home (1985), we adopt the 
map of maximum entropy, which can loosely be thought of 
as the smoothest map, or map containing the least infor- 
mation, with respect to a default map which may contain 
some a priori information about the surface intensity distri- 
bution across the donor star Generally, we assume a uni- 
form default map whose pixels are set to the average value 
in the reconstructed map, and the default map is updated 
after each iteration. An efficient algorithm for maximising 
entropy subject to a constraint has been implemented by 
Skilling & Bryan (1984) in the FORTRAN package mem- 
SYS. 

3 Roche Tomography highlights 

The earliest work using Roche tomography concentrated on 
reconstructing surface maps of CV donors using either sin- 
gle absorption or emission lines obtained using intermedi- 
ate resolution spectrographs. The maps clearly showed the 
varying impact of irradiation from the compact object and 
accretion regions on these stars. While interesting in their 
own right, these results have been covered in several papers 
and reviews already and, instead of repeating the findings 
here, we refer the reader to the work by Watson et al. (2003). 
Instead, we shall focus on some of the more recent results 
showing the starspot distribution on a number of C V donors. 

3.1 AE Aqr - the first image of starspots on a CV 
donor star 

Early, single line CV studies failed to reveal the presence of 
cool starspots on C V donors since they neither had the reso- 
lution nor the signal-to-noise to detect the small, faint signa- 
tures of starspots. (One should consider that most CVs have 
rriy >11, generally about 2-4 mags fainter than most stars 
that have been Doppler imaged). The first of these prob- 
lems can be solved by employing a high resolution echelle 



spectrograph. Simply moving to a larger aperture telescope, 
however, does not solve the signal-to-noise problem. 

Fortunately, the wide wavelength coverage of an echelle 
spectrograph allows lOOO's of stellar lines to be observed 
in one spectrum. One can then apply a technique known as 
Least Squares Deconvolution (LSD - see Donati et al. 1997) 
which effectively combines all the lines to form one 'aver- 
age' stellar absorption line of greatly increased signal-to- 
noise. Typically, this technique provides a multi-plex gain 
in signal-to-noise of ~30. To put this in perspective, using 
LSD on an 8-m class telescope is the equivalent of attempt- 
ing a single line study on a 40-m class telescope! 

Armed with this technique, we obtained two nights echelle] 
spectroscopy with the 4.2-m William Herschel Telescope 
(WHT) on La Palma on AE Aqr as part of a pilot study. This 
target was chosen as it is the brightest CV in the Northern 
sky, and the donor star contributes a large fraction (^ 70%) 
to the total system hght. AE Aqr also has a relatively long- 
orbital period of 9.87-hrs which meant that reasonably long 
exposures could be used before the donor's orbital motion 
started to smear out the fine profiles. Application of LSD to 
this dataset resulted in a multi-plex gain in signal-to-noise 
of ~26 over single-Une studies. 

The Roche tomogram of AE Aqr is published in Wat- 
son et al. (2006), and can be viewed as a movie on AN's 
webpages. Representative frames are displayed in Figure 1 
and show several dark starspot features. The most promi- 
nent of these is a large, high latitude spot most clearly vis- 
ible at orbital phase 0=0.25. This spot is centred on a lati- 
tude of ^^65°, in stark contrast to our Sun where few spots 
are seen above latitudes greater than ^35°. Such a high lati- 
tude feature on AE Aqr is, however, consistent with Doppler 
images of rapid rotators which often show the presence of 
polar spots. 

Another feature of note is the dark region surrounding 
the Li point, which we have interpreted as irradiation of 
the inner face by the white dwarf and/or accretion regions. 
Such irradiation zones in similar locations in CV donors 
have been mapped previously (e.g. Davey & Smith 1996; 
Watson et al. 2003). There does, however, appear to be a 
chain of spots extending from the polar regions down to the 
irradiated Li point, making longitudes of the star facing the 
white dwarf companion appreciably more spotted than else- 
where. Indeed, in the Doppler image of the pre-CV V471 
Tau, Hussain et al. (2006) also found that the side of the star 
facing the white dwarf was heavily spotted. This may well 
be due to tidal forces, which have been predicted to cause 
spots to emerge at preferred longitudes (e.g. Holzwarth & 
Schiissler 2003). 

In total, Watson et al. (2006) estimated that 18 per cent 
of AE Aqr's Northern hemisphere was covered with spots, 
providing the first conclusive evidence that these systems 
were magnetically active. This put the canonical theory of 
CV evolution (which required strong magnetic fields on the 
donors in order to drain angular momentum from the binary 
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Fig. 1 The Roche tomogram of AE Aqr. The panels (clockwise 
from top-left) show the cool donor star at orbital phases <l)=0 (infe- 
rior conjunction of the donor), <!)=0.25, (j) = 0.75, and finally 0=0.5 
(superior conjunction of the donor star). Dark greyscales indicate 
regions that are either covered in starspots or are irradiated. The 
system is plotted as the observer would view AE Aqr at an orbital 
inclination of i = 66° . 

and drive the systems to shorter orbital periods) on a more 
secure observational footing. 

3.2 BV Cen - starspots and slingshot prominences 

Based on the success of the AE Aqr pilot study we began a 
campaign on the 6.5-m Magellan Clay telescope situated in 
Chile to map other CV donors. The larger aperture, wider 
spectral coverage and higher throughput of the Magellan 
setup compared to the WHT allowed systems 2 or 3 mag- 
nitudes fainter then AE Aqr to be imaged. The second CV 
we observed was the ~14.7 hour period system BV Cen. 
A movie of BV Cen, reconstructed after applying LSD, is 
available on the AN website and some snapshots are shown 
in Figure 2. 

The maps of BV Cen and AE Aqr (Figure 1) make for 
interesting comparison. Again, BV Cen shows a high lat- 
itude spot at ~65°. Interestingly, the high latitude spot in 
both systems are displaced in the same direction towards the 
trailing hemisphere of the star (the side facing away from 
the orbital motion of the star). While we believe that the 
prominent dark region near the Li point is due to irradia- 
tion, we also see a chain of spots extending from the polar 
regions down to the Li point similar to that seen on AE 
Aqr. This provides some evidence that a mechanism is at 
work, perhaps tidal and/or Coriolis forces, that is forcing 
flux tubes to arise at these locations. 

The fact that regions near the Li point seem heavily 
spotted, combined with the fact that low latitude spots are 
seen near the Li point, has interesting implications for the 
accretion dynamics of these systems. Many accreting bina- 
ries show sudden dips in their lightcurves, and it has been 



Fig. 2 The same as Fig. 1, but for BV Cen and assuming an or- 
bital inclination of i = 53° . While B V Cen clearly has a larger spot 
coverage than AE Aqr, the similarities in the surface distribution 
of spots between BV Cen and AE Aqr is striking. 



suggested that this is caused by starspots moving across the 
Li point and quenching the mass transfer, leading to a dim- 
ming in the system's brightness. In their modelling of the 
mass transfer history of AM Her, Hessman, Gansicke & 
Mattel (2000) concluded that if starspots were to cause such 
dips then regions near the Li point would most likely be 
unusually heavily spotted, with a spot-filling factor around 
50%. Certainly, both AE Aqr and BV Cen appear to be more 
densely spotted on the hemisphere facing the white-dwarf. 
Indeed, we have calculated that, for B V Cen, the spot-filling 
factor reaches 40% at longitudes near the Li point, which 
seems to support the conclusions of Hessman et al. (2000). 

In addition to the presence of starspots in B V Cen, which| 
we calculated covered some 25% of the Northern hemi- 
sphere, we also detected a transient narrow emission feature 
at zero velocity. Such emission has been seen previously 
in other CVs (see e.g. Steeghs et al. 1996), also situated at 
zero velocity, and have been interpreted as 'slingshot promi- 
nences' from the donor star sitting near the centre-of-mass 
of the binary. The preference for prominences to be seen 
at this location may be a combination of two things. First, 
prominences erupting from near the Li point will be more 
subject to illumination from the accretion regions, making 
them more visible to us. Second, it appears from the Roche 
tomograms of BV Cen and AE Aqr that spots are more 
likely near the Li point, which may increase the probability 
of a prominence erupting in these active regions. Further- 
more, with such a strong concentration of activity near the 
Li point, magnetic fragmentation of the accretions stream 
may be responsible for the 'blobby' accretion seen in some 
CVs (e.g. Meintjes 2004; Meintjes & Jurua 2006). 
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Fig. 3 The same as Figures 1 and 2 but for V426 Oph, assuming 
an orbital inclination of 63° . 



3.3 V426 Oph - caught during a prolonged low-state 

Finally, we present a preliminary Roche tomogram of V426 



4 Future opportunities 

While we are now beginning to reveal in some detail the 
magnetic nature of CVs, a sustained campaign is required 
to investigate whether these systems show solar type activ- 
ity cycles, and to what extent such a cycle impacts on the 
observed behaviour of the binary. We are currently under- 
taking a monitoring campaign to determine whether AE Aqr 
displays any solar-like activity cycle. 

Another interesting question is whether the donor stars 
in interacting binaries display an appreciable level of dif- 
ferential rotation. Scharlemann (1982) suggests that tidal 
forces should weaken, but not suppress, differential rota- 
tion. The fact that we see very similar spot distributions on 
AE Aqr and BV Cen suggests that the spot distribution may 
essentially be fixed. In turn, this suggests that differential 
rotation on these systems may be weak, perhaps even sup- 
pressed. Given that differential rotation is a key ingredient 
in stellar dynamo models, amphfying and transforming ini- 
tial poloidal field into toroidal field, any suppression of dif- 
ferential rotation may provide a challenge to dynamo the- 
ory. Thus measurements of the differential rotation rate, via 
cross-correlation of latitude strips from images taken 1 or 



Oph, again as a movie with some representative frames shown|2 weeks apart, is essential. Such studies on binaries have a 



in Figure 3. With an orbital period of ~6.85 hrs this is, to 
our knowledge, the fastest rotating star on which starspots 
have been imaged to date, the previous record being 7.44 
hrs (RXJ1508.6 - Donati et al. 2000). 

UnUke the Roche tomograms of AE Aqr and BV Cen, 
V426 Oph does not show such a prominent high latitude 
spot, though there is still a hint of a weak polar feature. 



unique advantage over those conducted on single stars in 
that the sense of the differential rotation can be deduced 
since the actual mean rotation period is known from the or- 
bital period of the binary, allowing the corotation latitude of 
the donor star to be determined. 

Finally, it would be interesting to prove that starspots 
can affect the accretion dynamics of interacting binaries by 



again displaced shghtly towards the trailing hemisphere. V426fcuppressing mass transfer. An observation of a starspot cross-| 

Oph also shows spots appearing at lower latitudes. With a 
rotation rate almost 90 times solar, this is at odds with the 
models of Schussler et al. (1996) which predict that mag- 
netic flux tubes should only break the surface at high lati- 
tudes due to Coriolis forces and magnetic buoyancy. 

Although the reason for the lack of a high latitude spot 
is unknown, one possibility is that we caught it at an un- 
usual stage in its activity cycle. Interestingly, data from the 
AAVSO indicates that V426 Oph was undergoing some vin- 
usual behaviour at this time - being both dimmer during 
2005 (when this data was taken), and exhibiting far fewer 
outbursts compared to other years. Was V426 Oph in an un- 
usual state of activity, and is this responsible for its unusual 
behaviour during 2005? Further observations are required 
if any such hnk is to be made, or if we are to determine 
whether V426 Oph simply lacks a polar spot at all times. 

Other than that, the inner face of V426 Oph appears to 
be quite uniformly irradiated (giving the star a darker ap- 
pearance in the Roche tomograms at $=0.5). Within the ir- 
radiated zone, however, individual starspots are still visible 
near the Li point, and V426 Oph also appears to exhibit a 
less distinctive chain of spots extending from high latitudes 
down to the ii point. Given the preliminary nature of the 
analysis presented here, we are hesitant to draw any more 
conclusions at this time. 



ing the Li point, combined with a dip in the luminosity of 
the system, would be compelling evidence for such a sce- 
nario. It is quite comforting to see that, in the case of V426 
Oph, we can still identify starspots even if they lie within ir- 
radiated regions. In sunmiary, such work promises to reveal 
much about both the behaviour of binaries, and of the stellar 
dynamo under the most extreme conditions. 
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